COMMENT 



Lifshitz-point critical behaviour to 0(e 2 ) 

H. W. Diehl and M. Shpot 

Fachbereich Physik, Universitat Essen, D-45117 Essen, Federal Republic of Germany 

Abstract. We comment on a recent letter by L. C. de Albuquerque and M. M. 
Leite (J. Phys. A: Math. Gen. 34 (2001) L327-L332), in which results to second 
order ine = 4 — d + ^ were presented for the critical exponents z^L2, ?7L2 and 7L2 
of (i-dimensional systems at m-axial Lifshitz points. We point out that their results 
are at variance with ours. The discrepancy is due to their incorrect computation of 
momentum-space integrals. Their speculation that the field-theoretic renormalization 
group approach, if performed in position space, might give results different from when 
it is performed in momentum space is refuted. 
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In a recent letter |] de Albuquerque and Leite (AL) presented results to second order in 
e = 4 — d + y for the critical exponents 1^2, Vh2 and 7L2 of (i-dimensional systems with 
m-axial Lifshitz points. For the special case m = 1 of a uniaxial Lifshitz point, these 
results were previously given in a (so far apparently unpublished) preprint 0. The e 2 
terms AL found are at variance with ours |3|. [|] . 

As an explanation for these discrepancies AL suggest the following. Both we as well 
as AL employed a field-theoretic renormalization group approach based on dimensional 
regularization. To compute the residua of the ultraviolet poles at e = 0, we found 
it convenient to perform (part of) the calculation in position space. By contrast, AL 
worked entirely in momentum space. They speculate |I| "that calculations performed 
in momentum space and coordinate space are inequivalent, as far as the Lifshitz critical 
behaviour is concerned". 

This speculation is untenable and a serious misconception, a fact which should be 
obvious not only to readers with a background in field theory. The reason simply is: at 
each step of the calculation one can transform from momentum space to position space 
and vice versa. 

To become more specific, consider an A-point vertex function T^ N '(xi, . . . , x^) 
of a dimensionally regularized, translationally invariant, renormalizable Euclidean field 
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theory, such as the |0| 4 theory for an m-axial Lifshitz point considered both by us 
and AL. In momentum space the vertex functions have the form 

/ N 

f w (fli,..,fe) {^) d 5 J> 

\i=i 

where Y^ N > are conventional functions of N — 1 independent momenta, e.g., q 1; . . . , 9jv-i- 
The r^- 1 also depend on e (i.e., on d): they are meromorphic in e, having ultraviolet 
(uv) poles at e = 0. Since the issue is the e expansion, these are the only poles we 
have to consider; possible other uv poles at special values e > need not concern us 
here. Likewise, we do not have to worry about possible infrared poles one encounters 
in perturbation expansions about the Lifshitz point for a fixed space dimension d, nor 
embark on a discussion of related subtle questions like the appearance of perturbatively 
non-accessible mass-shifts meromorphic in e, and on how these problems are avoided in 
massive, fixed-<i renormalization schemes. 

The Fourier back-transforms of the functions define generalized functions 

(distributions), which depend on the N — 1 difference variables 
1 < j < N — 1. The same applies to each individual Feynman integral contributing to 

r (7V) . 

Field-theoretical renormalization group approaches such as the ones on which 
AL's work and ours are based hinge on the possibility of introducing a well-defined 
renormalized theory by absorbing the uv singularities of the theory's primitively 
divergent vertex functions in a consistent manner through counter-terms that are local 
in position space. In order that the renormalization procedure can be interpreted as a 
re-parametrization, these counter-terms must have the form of the (local) interactions 
appearing in the original Hamiltonian, except for a finite number of admissible additive 
ones. Well-known mathematical renormalization theorems 0, || ensure that this is the 
case, order by order in perturbation theory. 

Central to the proofs of such renormalization theorems is the observation that 
the primitive uv singularities have a local structure in position space. It is precisely 
this property that is crucial for the renormalizability of the theory. It ensures that 
the counter-terms, computed to a given order of perturbation theory, provide the 
subtractions for all divergent subintegrations of the Feynman graphs of the next higher 
order that are required to cancel all those uv singularities that do not have the form 
of local counter-terms. Such nonlocal uv singularities occur indeed: for instance, the 
graph ^-C^ has momentum-dependent pole terms ~ e _1 (involving logarithms of 
momenta). These are due to the divergent subintegral )CZX j they do not have the 
form of local counter-terms but cancel upon making the appropriate subtraction for this 
subgraph. (This subtraction is produced by part of the one-loop counter-term oc 4 ; 
see, e.g., Sec. 3.B of Ref. 0.) Zimmermann's forest formula || clarifies precisely which 
subtractions have to be made for each individual Feynman graph. The locality of the 
counter-terms manifests itself in the fact that in the final subtractions which must be 
made for superficially divergent graphs the graph is shrunk to a point. 
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What we have just explained has been known for decades and can be found in 
standard textbooks on field theory. It is true that many authors for computational 
reasons prefer the momentum representation when explaining the renormalization 
procedure. Therefore the significance of the uv singularities' local structure in position 
space may escape the reader's attention if not properly emphasized. However, a very 
clear exposition of the importance of this locality is given already in one of the earliest, 
classic textbooks on renormalization |J. 

The renormalization procedure can be performed equally well in momentum or 
position space. Utilizing dimensional regularization in conjunction with minimal 
subtraction of poles is advantageous in that one does not have to worry about how the 
regularization scheme and the conditions for fixing the counter-terms translate upon 
Fourier transformation: the scheme can be applied equally well in the momentum or 
position representation. The upshot of these considerations is that there is no way that 
AL's and our calculation can be both correct. 

The source of the discrepancies between AL's work and ours can be traced back 
to the different results they find for the required two-loop integrals. For example, our 
result for the integral hip, k) defined in Eq. (3) of AL's paper || reads 

j a {m)k A jcj,(m)p 2 



16m (m+ 2) 2 (8 - m) 



+ 0(6 U ), (1) 



with 



2l0+m Y(-) r °° 



where 



Mm) = / di; v m ~ l &(v- m, <f ) , (2) 

r(2-f)r(f) 2 Jo v > 

7 / 7 x 5 r~j — (3) 

(2n) d *- m J (2vr) m p 2 + k* y ' 

is the scaling function associated with the free critical propagator in position space (cf. 
Eq. (13) of Ref. ||), at the upper critical dimension d* = 4 + y. Here e is a unit 
d* — m vector, while v is an arbitrarily directed m- vector. The integral j a (m) is similar 
to j^(m), except that its integrand has an additional factor v 4 . 

From AL's Eqs. (11) and (18), we can infer their result for I$(p, 0); it reads 

4 AL )(p, 0) = -7T 4+ f V - + 0(6°) . (4) 

To see that this cannot be correct, one must merely consider the isotropic case 
m = d = 8 — e: for this, AL's result (f|) predicts a pole oc e -2 , even though the 
pole part oc p 2 must vanish because p has d — m = components. By contrast, our 
result does not violate this condition since j<f,(8) = 0. (See Sec. 4.5 and 4.4 of 

Ref. [|[ where we verified that our e-expansion results for general values of m reduce to 
known ones in both isotropic case m = d and m — 0, respectively.) 

AL realized the incorrectness of their findings for m — 8. Yet they seem to believe 
that the 'approximations' they made in their computation of i > 2 loop integrals do 
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not lead to erroneous results. Details of their approximations are described in Ref. 0. 
The crux of their method is 'to impose the constraint' ki = — 2k 2 on the momenta of 
the internal integral 

j, , . x f d d - m P2 d m k 2 

2(Pl+p ' l) ~J (Pi+kMpi+P2+P) 2 +(ki+k 2 n (5) 

of 

/Jd-m Amu 

p2 P l kf 1 / 2 (P 1 +P,fc 1 ). (6) 

This amounts to modifying the momentum term (fej + k 2 ) 4 of the last propagator 
in Eq. (|5|) to fcf- The error this introduces is given by the analogue of the integral 
(0) one obtains through replacement of I 2 (Pi +p,ki) by the corresponding difference 
SI 2 (., h) = I 2 (., fei) - J2O, 0), namely 



{pl + K)[ Pl + p 2 + pf + {k l + k 2 Y] 



k 4 + A + A; 2 ) fc! • fc 2 + 6 fc? fc 2 

[(Pi+P 2 +p) 2 + fc 2 4 ] • U 
Now the pole term oc p 2 /e of J3 we are concerned with corresponds to a logarithmic uv 
divergence ~ p 2 In A at the upper critical dimension (A = cut-off). In order for AL's 
approximation to be correct, SI 2 must have no contributions that vary as p 2 p^ 2 or p 2 k^ 4 
as p\ ~ k\ ~ A — > 00. As can be seen for instance by power counting, this condition is 
noi satisfied. (Readers preferring more mathematical scrutiny might want to compute 
Vp5I 2 and study its behaviour for large p± and hi.) Accordingly, AL's approximation 
is unjustified whenever m/0. The same kind of approximations are employed by AL 
for other i > 2 loop integrals. 

In closing, let us outline how the pole term oc p 2 /e of I3 given in Eq. (|6|) can 
be recovered via a momentum-space calculation. Using a Schwinger representation for 
each one of the three propagators in Eqs. and and performing the Gaussian 
integrations over p l and p 2 , we obtain 



d — m 
2 



/*OG POO POO 

hip, 0) = TT d ~ m / dx / dy / dz (xy + yz + zx) 
Jo io Jo 

Next, we make the variable transformations X = x/z, Y = y/z and K 12 = z 1//4 fc 12 , 
and take the derivative — d/dp 2 \ p=1 inside the integrals. The integration over z can 
now be performed; it produces the factor T(e) (1 + X" 1 + Y' 1 ) e = 1/e + O(e ). Upon 
transforming to the variables s = 1/X and t = 1/Y, one finds that 

-dl 3 (p,0) vr d *- m 



<9p 2 



/»oo /*oc 

ds dt(sty^ (l + s + t)f 
Jo Jo 



X 
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This is in conformity with Eqs. (|I])-(|3|). To see this, note that the integral is 
proportional to J d m t>$ 3 . In momentum space, this is a convolution of the form 
J k fca $fc 1 $fc 2 $fci+fc 2 - The Fourier transform can be read off from Eq. (14) of Ref. ||; 
it involves a modified Bessel function K u (k 2 ), for which we use the representation 

f°° i- 4 

k 2u K u (k 2 ) — 1 I dxx v - l e- x ~^ , (10) 
Jo 

with integration variables x, y and z. Employing the transformations s = x/z, t = y/z 
and Ki : 2 = 2 _1//4 fei i 2, we perform the integration over z. The result is the residuum of 
the pole (D- 

To summarize: AL's results are incorrect because of their unacceptable 
approximations made in computing Feynman diagrams. Their speculation that the 
field-theoretic RG approach might yield different results depending on whether it is 
performed in position or momentum space does not hold. 

We gratefully acknowledge the support by the Deutsche Forschungsgemeinschaft 
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